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We investigate electron transfer processes in donor-acceptor systems with a coupling of the elec- 
tronic degrees of freedom to a common bosonic bath. The model allows to study many-particle effects 
and the influence of the local Coulomb interaction U between electrons on donor and acceptor sites. 
Using the non-perturbative numerical renormalization group approach we find distinct differences 
between the electron transfer characteristics in the single- and two-particle subspaces. We calculate 
the critical electron-boson coupling a c as a function of U and show results for density- density cor- 
relation functions in the whole parameter space. The possibility of many-particle (bipolaronic) and 
Coulomb-assisted transfer is discussed. 



PACS numbers: 05.10.Cc (renormalization group methods), 71.27.+a (strongly correlated electron systems), 
82.39.Jn (charge transfer in biological systems) 
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Introduction - Electron transfer (ET) is a fundamental 
process in chemistry, biology and physics, for example 
in corrosion of metals, charge transfer in semiconduc- 
tors, enzymatic activities, cell metabolism, and photo- 
synthesis 0, 0- The characterization of ET processes 
in bio-molecules is an important step towards an under- 
standing of the biological function of many proteins and 
towards the possible construction of bio-molecular elec- 
tronic devices or biosensors. Theoretical investigations of 
ET processes typically start from a two-site model for the 
electronic degrees of freedom at the donor and acceptor 
sites which are coupled via a tunneling matrix element t. 
The correlated dynamics of electrons and vibronic modes 
(phonons) is essential for the ET characteristics. If the 
phonons are treated classically, one arrives at the Mar- 
cus theory Q . A quantum mechanical treatment of the 
phononic degrees of freedom J4J results in models related 
to the spin-boson model, in which the phonons are mod- 
eled by an infinite set of harmonic oscillators with a con- 
tinuous spectral density J (to) [UEJ. 

The following two limiting cases of ET processes are 
well understood. If the tunneling (or hopping) matrix 
element t is small we are in the limit of nonadiabatic ET. 
In the opposite limit, where t is large, the ET is adia- 
batic and the Born-Oppenheimer theorem holds. (The 
time of the electron moving from the donor to the accep- 
tor is too short for the vibronic modes to change their 
configuration.). Both limits are realized in ET processes 
in proteins but of particular interest are those param- 
eters which lie in the crossover regime. In this case, 
non-perturbative methods have to be applied (see, for 
example, Ref. 0). 

If the electronic part is treated in a one-particle pic- 
ture it reduces to two localized quantum states and can 
be modeled by a two-level system or spin, leading to a de- 
scription in terms of the spin-boson model. In many cases 
this may not be the adequate picture - many-particle ef- 



fects and electron-electron interactions have to be taken 
into account, for example if a more realistic modeling of 
the electronic degrees of freedom is required 0,0, if more 
than one electron is transfered simultaneously [lfj , or for 
a proper description of exciton transfer 

In this paper we propose a two-site electron-boson 
model to investigate many-particle effects and the role of 
electron-electron interactions. As sketched in Fig. fl] the 
model contains both a local Coulomb interaction U and 
the coupling to a common bosonic bath. Our approach 
goes beyond earlier work in which only the electronic 
degrees of freedom involved in the ET transfer were con- 
sidered (see Refs.HH[ll|) 

or the coupling to the bosonic 
bath was treated within the spin-boson model without 
considering many-particle effects and electronic correla- 
tions. 

We discuss competition or cooperation of vibronic and 
electronic effects (electron correlations) during the ET 
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FIG. 1: Schematic view of the two-site electron-boson model 
investigated in this paper. The tunneling matrix element of 
electrons between donor (D) and acceptor (A) sites is given 
by t; the parameter U denotes the local Coulomb interac- 
tion. Dissipation in the electron transfer process is due to the 
coupling of the electronic degrees of freedom to a common 
bosonic bath. 
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process by using the bosonic numerical renormalization 
group (NRG) method 0G3 

which is non-perturbative 
and can be used in the whole parameter regime between 
the adiabatic and nonadiabatic limits. We indeed find 
a significant difference between the ET characteristics in 
the single-particle subspace (which can be mapped onto 
the spin-boson model) and the two-particle subspace in 
which many-particle effects have to be considered. The 
coupling to the bosonic bath turns out to favor the for- 
mation of bipolarons, whereas a local Coulomb repulsion 
tends to delocalize the electrons. 

Model and Method - The two-site electron-boson model 
is given by the following Hamiltonian: 
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The operators c^J denote annihilation (creation) opera- 
tors for fermions with spin a on the donor (i = D) and 
acceptor (i = A) sites; riA/D is defined as n, = c L c i<j- 
The first three terms of the Hamiltonian eq. cor- 
respond to the two-site Hubbard model investigated in 
Ref. 8, with £j the on-site energies, t the hopping matrix 
element, and U the local Coulomb interaction for two 
electrons on either donor or acceptor sites. This part of 
the model can be easily extended to include non-local 
Coulomb correlations between electrons on donor and 
acceptor sites; for simplicity, we restrict the discussion 
here to local Coulomb terms only. The last two terms in 
eq. |QJ describe the free bosonic bath and the coupling 
between electrons and bosons, respectively. In the follow- 
ing we set 9a = 1 and </d = — 1 5 assuming a symmetric 
shift of the phonon displacements due to the electronic 
occupancy at donor and acceptor sites j2j. 

In analogy to the spin-boson model |a, ||| , the coupling 
of the electrons to the bath degrees of freedom is com- 
pletely specified by the bath spectral function 
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Here we assume an Ohmic bath, corresponding to the 
situation in many proteins 0], with J(u>) = 2irauj, < 
u> < u> c , with a the dimensionless coupling strength and 
lo c a cut-off which sets the energy scale in the following 
(w c = 1). 

The Hamiltonian eq. Q conserves both particle num- 
ber n e i and total spin of the electrons. The Hilbert space 
of the full model can therefore be divided into subspaces 
labeled by (Q,S Z ), with Q = n c \ — 2 (so that Q = cor- 
responds to half-filling) and S z the z-component of the 
total spin. 



Here we only consider the subspaces (Q,S Z ) = 
(-1,1/2) and (Q,S Z ) = (0,0). In the subspace (-1, 1/2), 
which we term the single-electron subspace, the model 
eq. m is equivalent to the spin-boson model 
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(3) 

with t = A/2, e A = e/2, and e D = -e/2. The two- 
electron subspace (Q,S Z ) = (0,0) is the main focus of 
the calculations presented in this paper; The electronic 
degrees of freedom in this subspace can be represented 
by the four-dimensional basis 



N) = {|U,o),|U>,II,T),|o,n>} , 



(4) 



with the notation \a, d) describing the occupation at the 
donor (d) and acceptor (a) sites. Consider now the 4 x 
4-matrix M c ^ =0 = (i\H ch \j) = 1,...,4) with the 
matrix elements taken only with respect to the electronic 
degrees of freedom. Introducing the notation 

n n 

we arrive at the matrix 
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The matrix defines the starting point for our nu- 
merical calculations. 

The technique we are using here, the bosonic NRG, has 
been described in detail in Refs. 1X3, |lJ3 in the context of 
the spin-boson model. The basic features of the bosonic 
NRG are as follows: (i) the logarithmic discretization of 
the bath spectral function J(u>) in intervals [A _n+1 , A~"], 
with n — 0, 1, . . . , oo and A > 1 the NRG discretization 
parameter (all the results shown in this paper have been 
calculated using A = 2); within each of these intervals 
only one bosonic degree of freedom is retained as a rep- 
resentative of the continuous set of degrees of freedom, 
(ii) The mapping of the resulting Hamiltonian onto a 
semi-infinite chain, (iii) The iterative diagonalization of 
the chain-Hamiltonian via successively adding one site to 
the chain. 

The bosonic NRG has been shown to give very accu- 
rate results for the spin-boson model 0, One of 
its strengths is the flexibility to handle a variety of mod- 
els involving the coupling of a small subsystem (here the 
electronic degrees of freedom at donor and acceptor sites) 
to a bosonic bath. The application to the two-electron 
subspace of the two-site electron-boson model is therefore 
straightforward. 
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FIG. 2: Phase diagrams of the two-site electron-boson model 
in the zero-bias case, e = 0; a) dependence of the critical 
dissipation strength a® =0 and Q ?=" 1 on the Coulomb inter- 
action U for t = 0.1 in the two-particle subspace (solid line) 
and in the single-particle subspace (dashed line) , respectively, 
b) dependence of a?=° on the hopping matrix element t for 
various values of U . 

Results - Let us first consider the zero-bias case £a = 
£d = for temperature T = . In the single-electron 
subspace, this gives e = for the corresponding spin- 
boson model which allows the observation of a quantum 
phase transition between a localized phase and a delo- 
calized phase at a critical a c (A) > 1 0, 0- For the 
value of t = 0.1 wc find « 1.2, naturally indepen- 

dent of U for the single-electron subspace (dashed line in 
Fig. Ek). Note that a^" 1 deviates from the exact value 
a|? =-1 (i — > 0) due to the finite t and the value of A as 
discussed in detail in Rcf. 

In the two-electron subspace we also observe a transi- 
tion between a localized and a delocalized phase (solid 
line in Fig. However, the dynamics of the two elec- 
trons is now correlated which leads to a [/-dependent 
a®= (U). Note that even for U = we find o^ =0 ([/ = 
0) 7^ a|? =_1 . This is because the coupling to the common 
bosonic bath generates an effective interaction between 
the two electrons [the energy of the total system is lower 
when both electrons occupy the same lattice site (D or 
A) as compared to the case of equal occupancy on D and 
A which does not result in a displacement of the oscilla- 
tors]. In other words, many-particle effects are present 
in the system even if the Coulomb interaction is assumed 
to be very small which is in close analogy to the physics 
of the Holstcin model In our model, the situation 

is more complicated; due to the coupling to a continuous 
bath, the effective interaction is always frequency depen- 
dent (there is no limit in which it could be reduced to 
a static one). Nevertheless, the net sign of the effective 
interaction is clearly negative, in agreement with the be- 
havior of the double occupancy shown in Fig. [3] which 
increases with increasing coupling to the bosonic bath. 




FIG. 3: Density-density correlation function C(uS) as a func- 
tion of the frequency uj for (a) different values of U and fixed 
a — 0.5 and (b) different values of a and fixed U — 1.0 (t = 0.1 
in both figures). In the insets of (a) and (b), the double oc- 
cupancy is displayed as a function of U and a, respectively. 



This also explains, to some degree, the shape and the 
asymmetry of the phase boundary a® =0 (U): a repulsive 
U > partly cancels the attractive interaction from the 
coupling to the bosonic bath so that increasing U favors 
delocalization of the electrons (for a similar feature in the 
Hubbard-Holstein model, see Ref . UtI) . 

An interesting consequence of this increase is that the 
lines aj? =_1 (£/) and a® =Q (U) cross at some value of 
U* w 2.4; for U > U* there exists a region in the phase 
diagram with a^ =_1 (C7) < a < a^ =0 (U) where the sys- 
tem is localized in the single-electron subspace but de- 
localized in the two-electron subspace. This can lead to 
interesting effects of the electron transfer dynamics: the 
addition of a second electron to the donor site drives 
the system from the localized to the delocalized phase so 
that either one of the electrons is enabled to hop to the 
acceptor site (such a process could be termed Coulomb- 
assisted transfer) or both electrons hop in a single event 
(bipolaronic transfer). 

The dependence of the critical a c on t in the two- 
electron subspace is shown in Fig. \Bp for various values 
of U . As expected, increasing t tends to delocalize the 
electrons, similar to the behavior in the single-electron 
subspace (see, for example, Fig. 8 in Ref. 0). 

To study the dynamics of the electron transfer process 
in the two-electron subspace we calculate the density- 
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C(t) dt 



density correlation function C{ui) 
with 

C(t) = ~ {[n D {t) - n A (t), n D (0) - n A (0)}+) . (7) 

In the single-electron subspace, the quantity corresponds 
to the well-studied spin-spin correlation function [fj ■ 

Figure |3| shows the NRG results for C(ui) for var- 
ious values of a and U with the temperature set to 
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FIG. 4: Density-density correlation function C(u>) for U = 
1.0, a — 0.5, t = 0.1, and different values of the bias e. 

T = 0. As expected from the corresponding behavior 
in the spin-boson model, the slope of C{lo) for small lu 
diverges upon approaching the transition from the de- 
localized side. This can be seen in Fig. |2t> where we 
show C(oj) for fixed U = 1.0 and increasing values of 
a < cvj? =0 « 0.55. Such an increasing slope in C(ui) cor- 
responds to a decreasing electron transfer rate k, with 
k — > for a — > aj? =0 . 

In the insets of Fig. [21 we show the double occupancy 
d = c?a = ("■AT n Aj.) = Jo C^dw (which is equal to 
d^> for zero bias). The double occupancy increases with 
increasing a; this is similar to the Holstein model, where 
the coupling to the phonons favors the formation of bipo- 
larons due to the attractive effective interaction. Increas- 
ing the value of U, on the other hand, reduces the double 
occupancy as shown in the inset of Fig. |2Ji. As shown in 
the main panel of Fig. [3Ji, increasing U also leads to an 
increase of the electron transfer rate [a decreasing slope 
inC(w)]. 

As in the spin-boson model, no quantum phase tran- 
sition is observed for any finite e where the system is al- 
ways in the delocalized phase The dependence of the 
density-density correlation function on the bias e is shown 
in Fig. 2J At small frequencies, we again find C(ui) a w 
with the slope decreasing with increasing e. Here one 
cannot deduce directly information on the transfer rate 
because of its non-monotonic behavior as a function of e 

a 

Conclusion - To summarize, we have shown that many- 
particle effects can lead to significant changes of the 
ET characteristics in donor-acceptor systems when more 
than one electron is present in the system. In the two-site 
electron-boson model proposed here, these many-particle 
effects are due to both the explicit inclusion of a local 
Coulomb interaction U and an effective interaction in- 
duced by the coupling to a common bosonic bath. We 
found that, depending on the model parameters, the ad- 
dition of a second electron can lead to a transition from 
localization to delocalization, or vice versa, see Fig. [21 
The results for the density-density correlation function 
and the double occupancy show that increasing the cou- 
pling to the bosonic bath favors the formation of bipo- 
larons, while increasing the Coulomb repulsion U has the 
opposite effect. 



The bosonic NRG is a very flexible tool and we are 
planning to investigate a number of extensions of the 
present model in the future. A generalization of the level 
structure allows, for example, the modeling of excitation 
transfer (see Ref. Inland the models proposed in Ref. H^) . 
Furthermore, the extension to three or more sites allows 
the investigation of multistep transfer processes and the 
influence of dissipation in the transport through (short) 
one-dimensional chains (see, for example, Ref. Il9r ). In 
general, the method presented here is a starting point 
for a more realistic description of electron and excitation 
transfer in those molecular systems where many-particle 
effects have to be considered. 
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